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Warm inflation scenarios are studied with the dissipative coefficient computed in the equilibrium 
approximation. Use is made of the analytical expressions available in the low temperature regime 
with focus on the possibility of achieving strong dissipation within this approximation. Two different 
types of models are examined: monomial or equivalently chaotic type potentials, and hybrid like 
models where the energy density during inflation is dominated by the false vacuum. In both cases 
dissipation is shown to typically increase during inflation and bring the system into the strong 
dissipative regime. Observational consequences are explored for the amplitude of the primordial 
spectrum and the spectral index, which translate into constraints on the number of fields mediating 
the dissipative mechanism, and the number of light degrees of freedom produced during inflation. 
This paper furthers the foundational development of warm inflation dynamics from first principles 
quantum field theory by calculating conservative lower bound estimates on dissipative effects during 
inflation using the well established thermal equilibrium approximation. This approximation does 
not completely represent the actual physical system and earlier work has shown relaxing both the 
equilibrium and low temperature constraints can substantially enlarge the warm inflation regime, 
but these improvements still need further theoretical development. 
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I. INTRODUCTION 



Cosmological observations are consistent with an early period of inflation in the evolution of the Universe, which 
among other things accounts for the observed flatness of the Universe and give rise to the primordial curvature 
perturbation. The quest now is for a realistic particle physics model of inflation, and a better understanding of 
the inflationary dynamics. The model has also to include the mechanism for the subsequent transition from an 
inflationary universe into a radiation dominated one. This requires the inflaton to couple to other particles, to allow 
for the conversion of the vacuum energy into radiation at the end of inflation during the reheating period. In addition, 
these interactions can lead to dissipative effects already relevant during inflation, affecting the inflationary predictions. 
These kind of scenarios are called warm inflation [lj, |2j, |3( . 

In the warm inflationary models, the dissipative term appears as an extra friction term in the evolution equation 
for the inflaton field cj>, and the corresponding source term for radiation pr\ 

4> + 3H<j> + r^ + V^ = 0, (1) 
p R + AH p R = T^ 2 , (2) 

where V<j, denotes the derivative of the potential with respect to the inflaton field. Warm inflation scenarios exhibit 
several features that are attractive for model building. For one, dissipation allows the inflaton mass during inflation 
to be much bigger than the Hubble scale 0, |{| , thus completely avoiding the "eta problem" [H 0, H, El IH3| i which is a 
generic problem in standard cold inflation scenarios in supergravity (sugra) theories. Another attractive model building 
feature is for monomial potentials, inflation occurs with the inflaton amplitude below the Planck scale, mp = 2 Ax 10 18 
Gev am. In contrast, for monomial potentials in the cold inflation case, usually called chaotic inflation scenarios 
[Tl|. the inflaton amplitude during inflation is larger than the Planck scale. This is a problem for model building, since 
in this case the infinite number of nonrenormalizable operator corrections, ~ X^^LiS™*^ n (W TO f )™ would become 
important and so have to be retained [To| . 

The challenge of realizing warm inflation is in understanding the dissipative dynamics from first principl es q uantum 
field theory. In terms of analytic approximations, much work has already been devoted to this problem |l2l ll'l IhH 
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E1E1 Inevitably warm inflation dynamics is nonequilibrium and a complete treatment of the statistical state during 
warm inflation is probably not amenable to simple analytic approximations, with more numerical based methods 
such as 0, 0] needed. Nevertheless, progress must be made systematically, initially understanding what can be 
from analytic and near analytic treatments. There are two known analytic approximations that would be relevant to 
apply to this problem. First is a quasiparticle approximation following Morikawa and Sasaki 19] . The other is the 
equilibrium approximation, which is based on the assumption that the statistical state of all fields always remains in 
thermal equilibrium pH l2*oj . The quasiparticle approximation has been developed for warm inflation in |l3l Il4l Il5| 
and several interesting model building possibilities have been found in 5J. The equilibrium approximation has been 
developed for warm inflation in |16( . In this paper we will explore some of the interesting warm inflation models that 
such an approximation yields. 

Testing the equilibrium approximation has an important general significance to the overall understanding of dis- 
sipative effects during inflation. In this approximation, the basic assumption is that the field system is minimally 
disturbed and so this approximation likely provides a lower bound on the degree of dissipative effects and particle 
production during inflation. Thus models that work under this approximation provide a minimal expectation on the 
overall robustness of warm inflation. Such solutions provide crucial existence proofs of the viability and consistency 
of warm inflation models with quantum field theory, which has been a basic question about these scenarios since they 
were first suggested 0,1^. This of course is of general significance in the development of inflationary dynamics, since 
an alternative way to state these results is regimes which would have been regarded unquestionably to be governed by 
standard cold inflation dynamics in fact are not, since dissipative effects are shown to be significant. For the better 
part of the existence of the inflation idea, the blanket assumption has always been the dynamics is cold inflationary. 
The observation was made much later in the development of the inflation subject p] that in fact this is not the 
unique situation and that dissipative effects can occur during inflation. Thus the results found in this paper provide 
a significant step in breaking these preset early notions about inflation dynamics and demonstrating that dissipation 
is a generic feature during inflationary expansion. 

In the high temperature regime, in earlier work a class of warm inflation models were found for the equilibrium 
approximation in |l2j| ; however one of the pronounced features of these models were the requirement of a very large 
number of fields ~ 10 4 , and so not attractive for most model building purposes (although see |22j for string motivated 
models). In this paper, we will show that in the low-temperature regime, many types of reasonable warm inflation 
models can be found, in particular for a moderate number of fields and sensible values of the parameters. 

The specific field theory models examined in this paper are the same as in |l3l Il4j , in which the inflaton field 
is coupled to a heavy bosonic field which in turn is coupled to light fields. An important feature of such coupling 
configurations is that even for large perturbative couplings, in supersymmetric (susy) models, the quantum corrections 
to the effective potential from these terms can be controlled enough to maintain an adequately flat inflaton potential 
|l4i Il5j. yet susy provides no cancellation of the time nonlocal terms, which are the dissipative terms, and so such 
terms can be quite significant. 

The basic picture of dissipative dynamics for this class of models is as the background inflaton field moves down 
the potential, it excites the heavy bosonic fields which in turn decays into light degrees of freedom EI Eld. The 
latter quickly thermalize and become radiation. Consistency of the approximations then demands the thermalization 
timescale to be smaller than the evolution timescale of the inflation field and the expansion timescale of the Universe. 
In addition, the condition T 3> H allows to neglect the expansion of the universe when computing T^,. Under these 
conditions, the dissipative coefficient T^, has been recently computed for a generic suspersymmetric inflationary model, 
in a thermal approximation 16] . It was shown that in the low T limit, with H <C T < m x , the dissipative coefficient 
behaves as T 3 /m x , m x being the mass of the mediating field. Still, dissipation can be large enough to dominate over 
the expansion rate and bring the inflaton into the strong dissipative regime with > H. 

In this paper we shall explore which kind of models of inflation can be brought into the strong dissipative regime 
within the low T approximation. We want to explore the viability of the scenario by checking whether we can satisfy 
the conditions > H and H < T < m x for at least 50-60 e-folds. Due to the T dependence on the dissipative 
coefficient, the ratio T/m x tends to increase during inflation, and at some point the low T approximation breaks 
down. The system would move then into the high T regime, where the dissipative coefficient goes linearly with T, 
and soon after that drops below H and inflation ends. The transition from the low to the high T regime may last 
some e-folds, but taking properly into account this period would require a more involved calculation of the dissipative 
coefficient beyond the analytical approximation. Therefore, as a first step in studying this kind of models, we only 
work within the low T approximation, assuming that inflation ends soon after the condition T < m x is violated. 

In section [n] we briefly review the basics of the warm inflationary dynamics in the strong dissipative regime, and 
applied this to some generic inflationary models, divided into two groups: (a) monomial potentials or chaotic inflation 
models, and (b) small field models where a constant term dominates the potential energy. In section lTTTl we summarize 
our findings and further comment on the consequences for model building. 
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II. DISSIPATION IN THE LOW T REGIME 



We consider the interactions given in the superpotential: 

W = g$X 2 + hXY 2 , (3) 

where $, X, and Y denote superfields, and </>, X; and V wm refer to their bosonic components. During inflation, the 
field y and its fermionic partner y remain massless, while the mediating field \ gets its mass from the interaction 
with the inflaton field <fi, with m x = V2g4>- Following Ref. the dissipative coefficient in the low T regime is well 
approximated by: 

T^0.64x s ¥^y^. (4) 



rriv I m 2 



Alternatively, we may consider a susy hybrid model of inflation, with the mass of the \ field given by m 2 . = 2g 2 (<fi 2 — (f> 2 ), 
with <p c being the critical value. In either case, the dissipative coefficient given in Eq. J3J does not depend on the 
coupling g in this regime, but only on h which can be quite large, h ~ O(l). The numerical coefficient in Eq. Q was 
computed taking X , Y to be singlet complex fields. But in principle, they may belong to larger representations of a 
Grand Unification Theory (GUT) group, as it is typically assumed in susy hybrid models. This will give rise to an 
extra factor of J\f = N x N% ecay in front of the dissipative coefficient, where N x is the multiplicity of the X superfield, 
and N decay counts the no. of decay channels available in X's decays. Taking m 2 , — 2g 2 <fi 2 , we have then: 

T 3 

~ CV»-^2 , (5) 

where C* = 0.16 x h A M. 

The conditions for slow-roll inflation are modified due to the extra friction term in Eq. Demanding <f> < 
(3H + T<p)(j) and <fi 2 < V, they are given now by: 

r ? = M I < 1 ' ( 6 ) 



(1 + r) 2 

- eH 
e - (i + r )2 

3H 2 (1+r) 3 



< 1, (7) 

< 1, (8) 



where rjn = mp(V^)/V, €h = m'j>(V < p/V) 2 /2 are the standard slow-roll parameters, and r — T ( p/(3H). In addition, 
once the source term for the radiation dominates in Eq. J2J, this would reduce to: 

AH PR ~ T^ 2 . (9) 

However, in order to have pr < AHpR and then Eq. ©, we need to impose instead the slow-roll conditions: 



(1 + r) 

eg 
(1 + r) 

3H 2 (1 +r) 



< 1, (10) 

< 1, (11) 

< 1- (12) 



For the kind of potentials considered in the following, we would have e# <C (V^ / 'phi) / '(3H 2 ) ~ r/H, and therefore the 
slow-roll conditions reduce mainly to: 

W ~ <1. (13) 



3iJ 2 (l + r) (1 + r) 



Once this is fulfilled, also the conditions r) < 1, e < 1 are satisfied. The above condition ensures that the radiation does 
not increase too fast during warm inflation, otherwise it will dominate too soon, not allowing inflation too proceed. 
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Therefore, in the strong dissipative regime with r > 1, the slow-roll evolution equations for the inflaton field and 
the radiation are given by: 

(14) 

PR - ahy; ■ (15) 

Alternatively, using pp = CrT , we have for the temperature of the thermal bath: 

where Cp — 7r 2 <7*/30, and is the effective number of light degrees of freedom. Once all the particles are in thermal 
equilibrium at a common T through rapid interactions, g* counts the number of relativistic degrees of freedom in the 
model. However, thermalization is not an instantaneous process, and for a particle specie kinetic equilibrium with the 
thermal plasma is only reached once its interaction rate becomes larger than the Hubble rate. Taking into account the 
thermalization rate would translate into a lower effective T than expected 23], or equivalently a lower value g*{T). 
Therefore, we may expect also g*(T) varying while in the dissipative regime, from say O(10) to the value for the 
Minimal Supersymetric Standard Model (MSSM) g* = 228.75 . 



A. Monomial potentials 



We first study the inflationary trajectory for a general inflaton monomial potential: 

-V 



V(4>) = V [-^\ , (17) 

with n > 0. Without enough dissipation, i.e., either for cold inflation with r = 0, or only weak dissipation with 
r < 1, these kind of models lead to inflation only for values of the inflaton field larger than the Planck mass mp. On 
the other hand, in the strong dissipative due to the larger friction term slow-roll conditions Eqs. © and Q can be 
fulfilled for values of the field well below Planck. That is, we can regard now Eq. 1)17(1 from the effective field theory 
point of view, with the potential well define below the cut-off scale mp; higher order term contributions suppressed 
by mp will be then negligible, without the need of fine-tuning the coefficients in front . 

During slow-roll, once the system enters in the strong dissipative regime, we can integrate Eq. I)14|l using Eqs. I(16|l 
and JSJl. The potential decreases from its initial value V^O) as: 

V ~'l-^0 7 . (18) 



V(0) V TV'(O) 

where N e is the no. of e-folds from the beginning of inflation, an pp(0) is the initial value for the radiation given by: 

1/7 

! , (19) 



p R (Q) ^ 9n 8 C% f 



V(0) \ 4 C\ \V(Q) 
Eq. I|18|l is only valid in the strong dissipative regime, r > 1 , with 

1/7 



^{swci) {W)J Uri {—) ' (20) 

which increases during inflation for n < 14. On the other hand, the values of r/n /r, T/H and T/<fi are given respectively 
by: 



(22) 



/ 81n 2 \ 1/7 / m% \ 2/7 (V(0p 2 7 



T 

H ~ \4CpCpJ \V{0)J V V 

T ( m P \( 81n 2 \ 1/7 fV(0)\ 3/U f V X 3 / 14 -V« 



" \V3<j>{Q)J \4C 4 C R J \m% ) \V(0) 



(23) 
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FIG. 1: Quartic potential: Allowed values for g* and Af = N x N^ ecay for having at least 50 e-folds of strong warm inflation: (a) 
enclosed region by the solid lines when V(0) 1/4 /m P = 0.3; (b) to the left of the dashed line when V(0) 1/4 /m P = 0.1. We have 
taken (f>(0) /mp = 1. 



Therefore, all the above ratios increase as inflation proceeds. This means that after some no. of e-folds, we will either 
have T/4> > 0.1, and the low T approximation will not longer hold, or T)h/t > 1 , and inflation will end shortly after 
when pb. > V. From (|21|l and 122(1 we have the relation: 

and for any power n > 1 we have that the radiation would not dominate the total energy density in the slow-roll 
regime with rju/r < 1. Once this ratio grows beyond one, so does the ratio Pr/V, and inflation ends when the 
radiation becomes larger than the inflaton potential. 

Without a specific model at hand, we can always consider Af = N x N% ecay and g* as free parameters 1 , and see for 
which values we can keep rin/r and T/<f> small enough for at least 50 e-folds or so (and T/H > 1). In turn, these 
values will depend on the initial value for V(0), and in the case of (f>/T also on (f>(0). For example, in Fig. we 
have plotted the allowed region in the plane g* — Af for a quartic potential taking V(O) 1 / 4 ~ 0.3mp and (f)(0) = mp; 
in this case, in order to satisfy all the constraints we require <?* < 100 but Af > 2300. Similar results are obtained for 
other powers of the potential. By lowering the value of the potential, it is easier to fulfill all conditions except that 
for the ratio rju/r- Keeping the latter below one gives the lower bound: 

AT > 8.4 x 10-VJ 4 (^) (n^(n - 1) + , (25) 

and the lower V(0) is, the larger Af has to be. For example, for n = 4, y(0)^ 4 /mp = 0.1, and — 10 we would 
need M > 2800 (dashed line in Fig. (JTJ), but getting to the no. of degrees of freedom for the MSSM, g„ = 228.75, 
would require Af > 29000, which looks rather large. In any case, curiously enough, due to the extra friction we can 
have inflation for values of the field below Planck, but the model prefers an initial value of the height of the potential 
only an order of magnitude or so below Planck. 

On the other hand, the amplitude of the primordial spectrum is also affected by the strong dissipative friction term 
and the presence of a thermal bath. Approximately, when T > H the fluctuations of the inflaton field are induced by 



We have also the value of the coupling h among the heavy field x an d the light degrees of freedom, which hereon we fix to its maximum 
allowed value, h = %/47r when giving bound on the value of Af. 
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the thermal fluctuations, instead of being vacuum fluctuations, with a spectrum proportional to the temperature of 
the thermal bath. In particular, when r>l we have for the spectrum of the inflaton fluctuations 0: 

<-(t) 1/4 ^< ™ 

with the amplitude of the primordial spectrum of the curvature perturbation given by: 




which in the case of the monomial potential is given by: 




In order to keep the amplitude of the primordial spectrum consistent with WMAP's value (24|, P^[ 2 ~ 5.5 x 10~ 5 , 
we would need a potential much smaller than O(10~ 12 m,p). But for such a value of the potential, we would need 
roughly Af ~ O(10 6 ) in order to get at least 50 e-folds in the strong dissipative regime, which seems rather large and 
unnatural 2 . Therefore, although in principle having strong dissipation during inflation is possible with a monomial 
potential, maintaining inflation for at least 50 e-folds requires rather large values of Af, and it tends to produce a too 
large amplitude for the primordial spectrum. 



B. Hybrid like models of inflation 

We consider now small field models of inflation, by adding a constant term to the monomial potential: 

n > (29) 

n = (30) 

Given that during inflation the potential is dominated by the constant term Vo, we can easily keep the value of the field 
below Planck in this class of models during slow-roll inflation. We regard this kind of potential as a generalization 
of a hybrid model 0, |2^|, where inflation ends once the inflaton field reaches the critical value, destabilizing the 
waterfall field coupled to it. Those interactions are not relevant to study the slow-roll dynamics, only to mark the 
end of inflation, and therefore we do not need to consider them in the inflationary potential Eqs. (|29|l . (|30|l . However, 
the same interactions between the inflaton and the waterfall field required by the hybrid mechanism will give rise to 
dissipation, and leads to Eq. (JSJ) in the low T regime. The case n — 2 would be the standard hybrid model [fil l2fij. 
with a mass term for the inflaton, whereas when n — we have the susy model with the logarithmic correction 
coming from the 1-loop effectivepotential |26l l27|. In supersymmetric hybrid models, one needs to worry about the 
so called "eta" problem [a S B S E3 > i- e -> the fact that generically sugra corrections give rise to scalar masses of the 
order of the Hubble parameter, including that of the inflaton, then forbidding slow- roll inflaton. Different solution 
to this problem exist in the literature, for example by combining specific forms of the superpotential and the Kahlcr 
potential [a 129- Nevertheless, typically although we can avoid the quadratic correction, i.e, a mass contribution, 
sugra corrections manifest as higher powers in the inflaton field 28]. In the case of strong warm inflation, the presence 
of the extra friction term alleviates the problem: slow-roll conditions are fulfilled also for inflaton masses in the range 
H < m<j, < T^. In addition, the values of the field being smaller than in standard cold inflation, the effect of higher 
order sugra corrections is also suppressed. 

The Hubble rate remains practically constant and given by H ~ Vq / '(y/3mp) during inflation, so neglecting its 
variation we can integrate the slow-roll equations in the strong dissipative regime, Eq. 114|) together with Eq. (TSJ, 



V(<t>) = V (l + 



M 



V{4>) = V 1 + In 



M 



2 Strictly speaking, the numbers and order of magnitude estimations have been obtained for the quartic monomial, but similar values are 
obtained for the quadratic and other powers. 
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and we obtain: 



( 4> 



T 



n ( Ma 2 C R 
2n 



1/7 



i/7 



-7/n 



T 



AC R C^ \H, 

(m \ ( <t> v- 2n/i 



and therefore: 



\T(o)j \m 



3/7 



VH 

r 

PR 

V 



3 \4C^C R 
(PR\ ( J>_ x 



H 



6n/7-2 



V V Jo V<Xo). 

where the subscript "0" denotes the initial value and we have defined: 

2 f (.^)(f) n . «^0, 



\ H 4 



71 = 0. 



(31) 
(32) 
(33) 

(34) 
(35) 
(36) 

(37) 



The amplitude of the primordial spectrum is given by Eq. (|27[1 . and the prediction for the spectral index ns is given 
by: 



n s - 1 =s 2 , . ? ~ 2- 



2V^ 



dlnB 



1/2 



dln/c 



3# 2 (1 +r) 



(38) 



evaluated at fco, corresponding to horizon exit (k — Ha) at say 50-60 e- folds before the end of inflation. Taking into 
account the T field dependence, derived from 4HpR ~ T^f/) 2 , we get: 



ns - 1 



f 



l + 7r 



3(3r-l) 6(l + 3r) m P 

-ejj - 3?7h H — —v2e H 



l + r 



l + r 



(39) 



where t/h and e# are the standard slow-roll parameters without dissipation. By taking again the derivative of Eq. 
(|39|l with respect to lnfc, i. e., with respect to the field, we obtain the running of the spectral index: 



, _ dn s 
Us ~ d\nk 



1 



1 - ns i 

r -i 

l + r 1 + 7r 

1 



3(3r - 1) £ , 



l + r 



H 



, 6(1 + 3r) m P - — 
3?ta H — 1— —V^h\ 



l + r 



(l + r)(l + 7r) 2 

with the derivatives of the slow-roll parameters and r given by: 



30e ff + 18j7h - 24[-fV2^} , 



>)'h 



e H = —- 



l + r 
2e„ 



(tih - 2e H ) ■ 



,mp 



H 



l + r 

mp\ ( vnp 

l + r 



'2e H + 2e H - m , 



(10e H + 8-^72^-6^), 



l + 7r 



(40) 

(41) 
(42) 

(43) 
(44) 
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where = 2mpV'" /V '. For the models considered in this section, we have that e# = (((j>/mp)r]H /(n — l)) 2 /2, and 
therefore typically e# <C r?jj, so keeping only the leading terms in Eq. (|39[l and (|4U|) . we end with 

ns-1 « ^( 7 -^+(±X , (45) 

and 

, ^ (m\ 2 ( <7-n) ( <j> V (14+10n-17n 2 ) \ 




(47) 

Notice that the spectral index if of the order of 0{r)H /r), whilst the running is of the order of 0((t]h /r) 2 ). Therefore, 
the same condition needed to have slow-roll in the strong dissipative regime will avoid having a too large spectral 
index. We will have a blue-tilted spectrum when n < 7, including the case of n = 0, i.e., the logarithmic potential. 
This is disfavoured by the data |24| when there is no running of the spectral index and no tensor contribution to the 
spectrum, otherwise the data is not conclusive. For example with non-negligible running, we have the allowed range 
0.97 < n,s < 1.21 0, with the running in the range —0.13 < dns/d\nk < 0.007. The more negative running, the 
more blue-tilted the spectrum can be, which would be the case for < n < 7 with n' s w — (ns — l) 2 n/ (3(7 — n)). 

Given that the field decreases during inflation, so does r/n/r, and also pr/V (or equivalently T/H) for any power 
n 7^ 0, being constant for the logarithmic potential n = 0. Therefore, the energy density in radiation will never 
dominate in this regime. On the other hand, <p/T diminishes for n < 4, but r only for n > 2. Therefore, for a 
logarithmic or quadratic potential once the system is brought into the strong dissipative regime stays there until the 
end. Indeed we may start in the weak dissipative regime, with T > H but r < 1, and it will evolve into r > 1. In the 
weak dissipative regime the inflaton evolution is given by the standard slow-roll equation, but the T behaves like: 



T C,,a A 



H 36C fl V# 
and then 



2(n-2) 

(48) 



3 \36CrJ \H 



(49) 



Therefore, for n — both T and r grow until T reaches a constant value when r > 1; and we may have the interesting 
situation of a transition from "cold" — > "weak warm" — > "strong warm" inflation. For n — 2, T remains constant until 
it diminishes in the strong dissipative regime, so that the parameter space divides into either cold or warm inflation, 
but the system can evolve from weak to strong dissipation. In the weak dissipative regime, quantum fluctuations of 
the inflaton field will also have a thermal origin, with an amplitude 0, 01 : 



pl/2 
TZ 



3H 3 



where in the second equality we have used Eq. Q48[l: the spectral index in given by: 

n s ~ 1 - 2 VH + 2e H , (51) 

Notice that even in the weak dissipative regime a logarithmic potential gives rise to blue-tilted spectrum, while for 
n > we would have a red-tilted spectrum, just the reverse than the standard cold predictions. 

For larger powers n > 2, we have the opposite behavior, and even if inflation starts in the strong dissipative regime 
it will evolve towards the weak and the cold regime. When this happens before the last 50-efolds of inflation, then 
dissipation becomes irrelevant. 

In the following we will briefly consider the cases n = 0, 2, 4 separately. The same than for the monomial potentials, 
we are interested in getting the constraints on the parameter values J\f and g*, which will depend on the power n, by 
demanding first having 50 e- folds of inflation with r > 1, and then getting the right primordial spectrum. 

Case n = 0: Hybrid logarithmic potential. For the potential Eq. I|30|l we have that the field value decreases 
exponentially as: 

4> ~ 4>(Q)exp[(r, H /r)N e ] , (52) 
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where the ratio rjn/r is approximately constant, and given by: 

m ( Ma 2 Cl 



r 



1/7 

(53) 



with a 2 — f3Vo/H 4 . The temperature as mentioned before is approximately constant and given by Eq. JSU with 
n = 0, while the ratios 4>/T and r decreases/increases respectively like: 

j, - I ~7j- ) ( a 4 ) exp[{ m /r)N e ] , (54) 
1 / H V fa 4 C 4 J 3 \ 3/7 

3 ( m ) ^£-J e Xp[ -(mMN e] (55) 

where we have used Eqs. (|33[1 . (|34[1 together with (|52[) . The lower bound on A/" is obtained in the limiting case for 
slow-roll warm inflation, when T/H ~ 1, r(0) ~ 1 and ((/)/T)n — 10, where the subscript "N" means N e e-folds after, 
which gives: 

TV ~ 11.87exp[~0.23N e gl /2 /TV 1/2 ] . (56) 

For example, with ~ 10, 7V e = 50 we have the lower bound: TV ~ 180, <j)(0)/H ~ 150, a 2 ~ 244, f^/r ~ -0.054; 
with 5* ~ 228.75, TV e = 50 we have: TV ~ 1350, <j>(0)/H ~ 1.2 x 10 3 , a 2 ~ 5.2 x 10 3 , and W ^ -0.1. Those values 
of TV = N x N% e are quite in the range of a realistic model; for example with x m the 126 or 351 of SO(10) {E e ), 
and Nj ecay w O(10), one can expect having TV in the range of a few thousands. However, for such values we always 
get a too large amplitude of the primordial spectrum, Eq. (|27|l . which using Eqs. (|53|l and l|55l) can be written as: 



>v» „ 0.4 r 1 



3/2 

Km J 1% 



^ -0.4 -- . (57) 



Therefore, in order to match the amplitude of the primordial spectrum with WMAP's value we need a larger initial 
value of the field 0(0), but then the value of Cp (TV) has to be larger in order to stay within the low T approximation, 
with 4>/T > 10; at the same time we need to increase a 2 in order to keep r(0) > 1. Satisfying WMAP's constraints 
requires then TV > 5 x 10 5 , rather large, with <j>{Q)/H > 2 x 10 6 and a 2 > 8 x 10 u . 

Having inflation in the weak warm regime, and a transition from weak to strong dissipation at the end, might in 
principle help in fixing the amplitude of the spectrum to lower values. We still need to impose that (a) T/H > 1, 
Eq. 148(1. (b) we get enough inflation, i.e., N e sa 50. These translates into TV > 0. 05 g */!]]], with rjn < l/(2N e ), and 
therefore we have the lower bound TV > 0.2g*N 2 . For example for 5* = 288.75 and N e ~ 50 we have TV > 1.2 x 10 5 , 
which again is rather large. Going from weak to strong dissipation earlier and having only 10 e-folds in the weak 
regime in principle is a possibility for smaller values of TV, but then again we cannot keep the low T approximation 
and 4>/T > 10 for the remaining e-folds of inflation. However one has to bear in mind that we are not properly taking 
into account how inflation ends, which requires going beyond the analytical approximations. The transition from the 
low to the high T is not necessarily instantaneous as we have implicitly assumed in this paper, and may last another 
few e-folds. Counting all together can bring again the values of the parameters into realistic values. 

Case n = 2: Hybrid Quadratic potential. The main restriction for this model comes from keeping the ratios 
T/H and <f>/T within the range of the low T approximations. During slow-roll the dissipative ratio r increases, so its 
initial value r$ gives the lower bound on this parameter. On the other hand, we still require {t)h/t)o = a 2 /(3^o) < 1 
in order to have slow-roll warm inflation. By using Eqs. I|31() . and (|34|l . we can rewrite Eqs. (|32|l . I|33() . in terms of ro 
and t)h / tq like: 

f^Ve)" 1 ) 2 , (58) 

3r (l + ^iV e ) 3 ) V2 , (59) 

Thus, unless Cr is too small, the condition <j)/T > 10 will be fulfilled if we have T > H. This in turn translates into 
a lower bound for TV: 

/V > 0.052$. + ^] . (60) 
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FIG. 2: Hybrid quadratic potential: Spectral index depending on the no. of e-folds to the end of inflation, for different values 
of N. We have taken: g* = 230, 4>(0)/m P = 0.21, r/ H = 3, and V 1/4 fmp = 3 x 10" 4 . 

From the approximated expression for the spectral index, Eq. (|45|) . if we want to keep ns within the observable 
range, we require nn/ro < 0.093, which for N e ~ 50 gives the lower bound Af > 32.5g». Again, having slow-roll warm 
inflation for example with ~ 228.75 needs Af > 7500, but for ~ 10 it only requires Af > 325. As an example, in 
fig. J2J we have plotted the predicted spectral index depending on the no. of e-folds left to the end of inflation, for 
Af = N x N% ecay — 10000, 20000, 30000 and g* = 228.75. The corresponding spectral index of the primordial spectrum 
would be that at around 50-55 e-folds, which is always ns < 1.2. The value of the running can be obtained from Eq. 
(pEtJj) . and it is given respectively by n' s ~ -2.5 x 10~ 3 , -8.5 x 10~ 4 , -4.7 x 10~ 4 . 
We can estimate the amplitude of the primordial spectrum, Eq. (|27|l in terms of the parameters r$ and rjjj/ro: 



Having a not too large amplitude for the primordial spectrum prefers values of tjh /tq as large as possible, but not too 
large values of ro. For values of Af, g* within the range of Eq. (|60[) . the amplitude remains below say 3 10 -4 for values 
of r$ of the order of O(10). Therefore in this kind of models parameter values can be found giving rise to the right 
order of magnitude for the primordial spectrum in the strong dissipative regime, but the stronger constraint comes 
from getting a not too blue-tilted spectrum. 

Case n > 4: Hybrid quartic and higher powers. In this case we have that T/H and r, Eqs. (|32fl and l|34|l . 

both decrease during inflation, while on the contrary <j)/T increases. Therefore, it is not difficult to remain within the 
range of the low T approximation, but on the other hand strong dissipation with r > 1 will last only a few e-folds. 
With n = 4, having 50 e-folds in the strong regime requires for example Af — N x N% ecay > 10 3 for ~ 10, and 
Af > 10 4 for g» ~ 228.75. In addition if we want to get the right amplitude for the spectrum Eq. I|27|) . and spectral 
index, this value increases by one order of magnitude, as we need to adjust the values of the quartic coupling a 2 = A 
to rather small values, and then the initial value of the field to larger values to have </>(0)/T > 10. Numbers do not 
change much if we demand 10 or 50 e-folds of inflation in the strong dissipative regime. 

Comments: Ending inflation 



3 Given the uncertainties in the analytical estimations, and in Eq. 1271 . it would not make sense to impose the specific WMAP value 
as a constraint, but an order of magnitude estimation would be sufficient. Furthermore, we have checked numerically that when the 
primordial spectrum is originated with ro not much larger than unity, the value of the amplitude tends to be lower than that given in 
Eq. 07J. 
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We have seen that depending on the values of the field (4>(0)) and couplings (Vo, 0) in the inflationary potential 
with n = 0,2, we may have either weak or strong dissipative regime during N e e-folds, but for rather large values of 
the parameter Af = N x N^ ecay . However, those large values are not required for having inflation in the weak/strong 
dissipative regime, but in order to match the predicted values of the spectrum with the cosmological observations, 
including having enough inflation in the low T regime. With lower values of the multiplicity parameter, in the range of 
a few thousands or so, dissipative effects are relevant for the inflationary dynamics, and indeed the tendency is to bring 
the system into the strong dissipative regime, first in the low T regime, and later into the high T regime, where the 
analytical approximations break down. From this point of view, cosmological observations rule out large regions of the 
parameter space in this kind of models. Nevertheless, a quantitative statement on the values of the parameters needs 
in turn to properly take into account the dynamics at the end of inflation. There are some important corrections 
that have to be included before ruling out any of the models presented here because they do not fit cosmological 
observations. First of all, as already mentioned, how the system interpolates between the low and the high T regime, 
and for how long (how many c-folds) this period lasts. In addition, in hybrid models as the inflaton evolves towards 
the critical value <j> c the mass of the waterfall field decreases accordingly, approaching the tachyonic instability, with 
m x = 2g 2 (<fr 2 — <fr 2 ) . The dissipative coefficient Eq. Q depends on the ratio of T 3 /m x , and thus close to the critical 
value it will get a rather large enhancement factor, with 



Although it may look at a first glance that this enhancement of T^, and therefore r, will bring the system faster into 
the high T regime, the increase of the extra friction in the inflaton evolution tends to slow down the field, and as a 
consequence also the evolution of the ratios T/H and <p/T. This again provides some extra e-folds of inflation in the 
warm dissipative regime. Finally, we will have to take into account thermal corrections in the scalar masses, mainly 
m x , which can be non negligible towards the end of inflation in the high T regime. Nevertheless, in order to take into 
account these effects one has to resort to numerical calculation, which are beyond the scope of this paper [30l |. The 
results presented here can be seen as a kind of worst case scenario, and the bounds on the parameters will be relaxed 
once we give up the requirement of having the full 50-60 last e-folds in the low T range. 



This paper has made an important step in understanding dissipative effects and particle production during the 
inflationary expansion period. These results are of general significance in understanding inflationary dynamics, since 
they break preconceived ideas set early in the development of inflation that the dynamics during inflation has negligible 
dissipative and particle production effects. The results in this paper show that picture does not express the general 
case and that dissipative effects become important in very generic models. Moreover as discussed in this paper and 
elsewhere 0, 0, H, 0, El El |Ml EI2 j the consequences of dissipation during inflation are nontrivial with respect to 
observable signatures and model building prospects. 

The technical problem of determining dissipative effects during inflation should not be underestimated. To appreci- 
ate this point, it is useful to compare this warm inflation problem to the standard cold inflation problem. In the latter, 
the particle production phase is pictured as entirely separated from the inflationary expansion phase; this problem 
has been intensely studied for well over two decades and is still not fully understood. The warm inflation problem 
is technically much more difficult, since particle production and inflationary expansion are meant to be occurring 
concurrently and this problem has been examined for a much shorter amount of time. 

The basic problem of determining dissipative effects during inflation can be posed as follows. A background inflaton 
field is evolving slowly in time as the Universe inflates. This field is coupled to other fields, thus in general dissipative 
and particle production effects can be expected. This is a nonequilibrium situation, with the central problem being 
to determine what is the statistical state of the Universe. To address this problem from first principles quantum 
field theory has required building up the necessary knowledge through an interplay of investigating nonequilbrium 
approximations and then testing their relevance to the actual inflation problem. Ideally one could imagine trying 
to solve this problem through numerical calculations and computer simulations, but the problem is too complex to 
immediately do a very general treatment this way. Ultimately the goal is to achieve some sort of reliable treatment 
of the problem by these methods. However before that can be achieved, a clearer idea of the approximate state of the 
system is needed. This is the direction that has been developing in the past few years The one baseline one has 
to go on is the assumption the statistical state remains close to thermal equilibrium throughout evolution. There are 
obvious deviations from equilibrium due to the evolution of the background inflaton field and due to production of 
light particles and their rescattering with the background field. Nevertheless, thermal equilibrium is a good limiting 
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case to examine since reliable and unambiguous calculations can be made of dissipation and particle production. 
However the assumption of equilibrium is a conservative one, from which one anticipates dissipative effects will be 
minimal. As such this approximation is vitally important in determining a minimal level of warm inflation. 

With these considerations in mind, the results in this paper are all the more encouraging. We have shown that 
all the simple inflation models, those with monomial and hybrid potentials, have warm inflation regimes. Moreover 
the two distinguishing model building features of warm inflation, inflaton mass bigger than the Hubble parameter so 
complete avoidance to the "eta problem" and field amplitudes below the Planck scale even in monomial potentials, 
have both been verified explicitly in models. Although improvements on the equilibrium approximation to obtain the 
correct nonequilibrium state will quantitatively change the results from those in this paper, these qualitative features 
will persist and so have been firmly established as realizable from realistic quantum field theory models. 

The quantitative changes from improving on the equilibrium approximation should lead to a decrease in the total 
number of fields and so to even more simpler models that yield warm inflation. In this paper, within the equilibrium 
approximation, we found models realizing warm inflation requiring a minimum of A^ x w O(10 3 /N% ecay ) fields. This 
is within the realm of realistic model building with particle physics GUT models such as 5*0(10) or Eg, with the 
X field in the 210 or 351, and a decay factor Ndecay ~ 0(3 — 4). To gain an idea of how much these models can 
be improved from the nonequilibrium treatment can be g ained by examining the results already studied for warm 
inflation using the quasiparticle approximation 0, Hj, [Hj • This approximation has been studied in the context 
of quantum field theory and cosmology applications since the work of Morikawa and Sasaki in the mid 1980s |19| . 
It is motivated by generic features known of many-body systems in condensed matter physics. No derivation of this 
approximation has been made in quantum field theory, but there is some suggestive numerical evidence in support 
of it |33| . One of the next steps in the development of warm inflation dynamics is to understand by how much the 
equilibrium approximation deviates and in particular how closely it tends to the quasiparticle approximation. 

In summary, this paper has accomplished two things. On the side of principle, due to the thermal equilibrium 
approximation underlying all the results, this paper has calculated minimal expectations of radiation production and 
dissipative effects during inflation. On the side of application, this paper has identified and developed a reliable 
methodology for performing warm inflation calculations. Although the results found here are only on the periph- 
ery of usefulness for model building, the methodology developed in this paper gives a foothold from which further 
improvements can be made. Aside from the major modifications already mentioned of expanding the equilibrium 
approximation to a more accurate nonequilibrium treatment, there are also several improves within the equilibrium 
approximation that can be made, and would lower the number of fields necessary, thus increase the scope of model 
building prospects. One modification is to move away from the strict low temperature regime into the intermediate 
temperature regime. This is technically much more difficult since now all finite temperature effective potential cor- 
rections will have to be accounted for. However the results of this paper already anticipate that the relaxation of this 
condition will increase dissipative effects, lower the total number of fields needed and help resolve issues of exiting 
the inflation epoch. A second improvement is treating higher order calculations of the dissipation coefficient. This 
dissipative coefficient is closely related to the shear viscosity, as observed in , for which resummation calculations 
|34| have shown provide as much as a factor four enhancement. One anticipates similar enhancements could occur 
also for the dissipative coefficients. 
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